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The characteristics and the scaling laws of isolated spherical blast waves have been 
briefly reviewed. Both self-similar solutions and numerical solutions of isolated blast waves 
are discussed. Blast profiles in the near-field (strong shock region) and the far-field (weak 
shock region) are examined. Particular attention is directed at the blast overpressure and 
shock propagating speed. Consideration is also given to the interaction of spherical blast 
waves. Test data for the propagation and interaction of spherical blast waves emanating 
from explosives placed in the vicinity of a solid propellant stack are presented. These data 
are discussed with regard to the scaling laws concerning the decay of blast overpressure. 


Nomenclature 


c 0 = ambient sound speed, ft/s 

d = distance from explosion, ft 

e = internal energy per unit mass 

E - total blast energy, lbf-ft 

f = frequency, Hz 

M = Mach number 

M s - shock Mach number 

N = number of coalesced waves 

p ~ static pressure, psia 

p 0 = ambient pressure, psia 

r - radial coordinate (distance from the source) 

R = shock radius, ft 

t = time, s 

U - fluid or particle velocity, ft/s 

U = shock velocity, ft/s 

W T = TNT equivalent weight of explosive, lb 
Z = Brode’s parameter, d / W r 1/3 (lb ,/3 ) 

Greek Symbols 

a =0,1,2 for plane, cylindrical, and spherical blast waves 

y - specific heat ratio 

X = d /(E / p 0 )' 13 
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p 

= 

density, lbm/ft 3 

4 


r/R(t) 


= 

c 0 2 //? 2 =l/M 2 

tip 

= 

shock overpressure, psi 

Subscripts 

0 

= 

upstream of shock (stationary) 

1 


downstream of shock 9stationary) 

C P 

= 

pressure coefficient 

Cx 

= 

force coefficient in the x direction 

Cy 

= 

force coefficient in the y direction 

c 

= 

chord 

dt 

= 

time step 

Fx 

= 

X component of the resultant pressure force acting on the vehicle 

Fy 

= 

Y component of the resultant pressure force acting on the vehicle 

f’g 


generic functions 

h 

= 

height 

i 

= 

time index during navigation 

j 

= 

waypoint index 

K 

= 

trailing-edge (TE) nondimensional angular deflection rate 


I. Introduction 

T he study of the propagation of blast waves in air is of great technical interest in many industrial applications- 
chemical, nuclear and aerospace. Blast waves are generated by point explosions (nuclear explosions and 
detonation of solid explosives, solid and liquid propellant rocket motors), high pressure gas containers (chemical 
explosions) and laser beam focusing in the ambient Blast wave problems also arise in astrophysics, hypersonic 
aerodynamics and hypervelocity impact. An understanding of the properties of the blast waves both in the near-field 
and the far-field is useful with regard to the characteristics such as blast overpressure, blast speed, and impulse. 

Interesting problems also arise in the propagation and interaction (coalescence) of blast waves emanating from 
the detonation of closely spaced explosive charges (sequentially detonating, spatially distributed ammunition 
stacks). Regions of significant overpressure associated with the coalescence of blast waves from distributed 
ammunition stacks may be less extensive than regulatory requirements (Starkenberg and Benjamin 1 ). 

Many theoretical and experimental studies were reported by several investigators on the planar, cylindrical and 
spherical blast waves since the classical work of Taylor 2 on strong (intense) blast waves from nuclear explosions. 
The decay of blast wave amplitude (pressure, density or other physical properties behind the shock) has been studied 
less often (Diaci & Mozina . This information will be useful in the possible detection of the source of the blast 
waves from far-field measurements. A proper understanding of the blast wave analysis and scaling is important in 
the interpretation of the test data. 

This report summarizes the blast wave scaling laws and the characteristics of isolated spherical blast waves. 
Test data for the propagation and coalescence of blast waves from closely spaced multiple charges and solid 
propellant are examined and analyzed. Procedure for Paper Submission 

II. Isolated Blast Wave Propagation 

When disturbances of finite amplitude are propagated in perfect gases (destitute of viscosity or heat conductivity), 
discontinuities in pressure, velocity and temperature of the medium may occur. These are called shock waves (or 
shocks) or blast waves (intense shock waves). Shock waves and the flowfield become planar, cylindrically 
symmetric, or spherically symmetric, respectively if the energy source is in a plane, along a line, or at a point in 
space. The appearance of shocks is a consequence of nonlinear character of the equations governing the propagation 
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of finite disturbances that is in part responsible for the formation of discontinuities. This nonlinear character is 
evident from the temporal evolution of the pressure wave (Fig. 1). 

Taylor 5,6 , in his classical blast wave theory, formulated the equations of motion in Eulerian form for the 
spherical blast waves. He assumed instantaneous release of energy is released in an infinitely small region (Fig. 2a), 
that the ambient pressure is negligible compared to the shock overpressure (i.e. strong or intense shock wave), and 
that the variation of transport properties are unimportant. The resulting similarity solution is valid for very short 
times after the blast occurs. 

Since the pioneering work of Taylor 6 , several treatises and reviews have been published, notably those of 
Sakurai , Sedov , Openheim et al. , Kinney and Graham . Sedov ’ and von Neumann obtained an exact closed 
form analytical solution to the blast flowfield problem considered by Taylor (the sol-called Taylor problem). This 
solution is valid for very short times after the blast (the Taylor problem). At later times the strong shock assumption 
breaks down, and numerical methods are required for a solution (Director and Dabora 13 ). 

Sakurai 14,15 presented a perturbation solution to the general blast wave problem. An exact numerical solution of 
the blast wave problem covering strong and weak shock regions is provide by Goldstine and von Neumann 16 . Back 
and Lee 17 proposed an improved analytical solution to this problem. 

Brode 18,19 numerically solved the Lagrangian equations in finite difference form with central differencing 
scheme. Both point sources and spheres of initially high pressure gas (Fig. 2b) and charges of TNT are considered. 


A. Governing Equations 
1 . Eulerian Formulation 

The equations of motion (continuity, momentum and energy) for the flow of 1-D, inviscid, compressible gas in 
Eulerian coordinates are (Director and Dabora 13 ): 



where r and t refer to space and time coordinates (independent variables), and e the internal energy per unit mass. 
The quantity a =0, 1, and 2 refers to planar, cylindrical, and spherical system respectively. The equation of state for 
a thermally perfect gas, supplementing the above equations, is expressed by 

p = pe{y-\) ( 2 ) 

where y denotes the isentropic exponent, and is assumed constant. 

2. Lagrangian Formulation 

The mass contained within any region of space adjacent to the origin is 


m = F 0 \pr a dr 

o 


( 3 ) 
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where F 0 =l, 2 7T , or 47T for plane, cylindrical or spherical system respectively. The conservation equations are 
transformed to Lagrangian form with the aid of the above equation as: 


|-(p r *)+F 0 pV ! ‘' 

dt dm 

d{e + u 2 / 2) i u d 

' * ft “I 


dm 
= 0 



(4a) 

(4b) 

(4c) 


3. Rankine-Hugonoit Relations 

Behind the shock, a strong blast wind is initiated, which moves in the direction of the shock. The particle 
velocity acquired from the shock waves is obtained from the Rankine-Hugonoit relations. The Rankine-Hugonoit 

relations (boundary conditions) across a shock of finite strength at r = /?(f)are (Shapiro 20 , p. 123; Liepmann & 
Roshko 21 , p. 64), see Fig. 3 : 
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Here the shock velocity U (relative to the undisturbed gas; i.e. relative to a fixed coordinate system) 

U = dR/dt (6) 

with R denoting the shock radius (or shock trajectory). The speed of sound in the undisturbed medium is defined as 

c 0 = 4wJ~Po (7) 


The shock Mach number M S (=U / c 0 =u 0 /c 0 ) and the particle velocity w^(near the shock) may be 
expressed in terms of the pressure ratio as 
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(8a) 


(8b) 


The quantity P = (p l /p 0 —\) represents the strength of the shock, and Ap = p x — p 0 denotes the 
overpressure. 

Strong Shock Conditions 


The strong shock conditions ( p x »p 0 , U»c 0 ) simplify to 
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2 


( 9 ) 


Pi _ r + 1 A _ 2/ U 2 Uj_ 

Po r - !’ /> 0 r + i« 0 2 ’ t 7 r + i 

In terms of the pressure ratio, the particle velocity may be expressed as 


2 Pi 

c 0 \r(r + i)poi 


( 10 ) 


B. Point Explosions 
1 . Strong Shock Solutions 

Self-similar solutions for the case of very strong shocks (short times after the explosion) were developed by 
Taylor 5,6 for spherical blast waves in his classical blast wave theory. These solutions are based on the form of the 
shock radius, which yield similarity solution. According to Taylor 2 , a finite but large amount of energy is suddenly 
released by nuclear fission in an infinitely concentrated form. This leads to the formation of a shock wave according 
to 


*(<)= S(r) P rE"V“ - SbiSEIpJ" ( 11 , 

which leads to a similarity solution. Here the quantity S is a dimensionless constant which is dependent on the 
specific heat ratio y , and is determined by solving the equations of motion. 

The governing partial differential equations (of continuity, momentum and energy) for a spherical wave (eqs. 
la-c in Eulerian form) are thereby reduced to an ordinary differential equation ODE, involving a single independent 
similarity variable defined by 


Z = r!R(t) 


( 12 ) 


Rankine-Hugonoit relations for strong shock (eq. 9) are utilized. The total energy (kinetic energy and the heat 
energy) behind the shock is expressed by 


= 4,J 


— pu 

2 H 


r 2 dr + 4 n 


Y - 1 


dr 


(13) 


which is assumed constant. The constancy of energy in each x — sphere implies that the energy flowing into the 
sphere with the new material that enters it exactly balances the work which its original surface does by expanding it 
against the surrounding pressure. The boundary conditions include the condition that the particle velocity at the 
center of explosion must be zero. 

Taylor 5,6 numerically solved the ODE. Taylor 2 solved the system for y - 1 .4 numerically, starting from the 
condition at the shock 7/ = 1 . An exact closed form analytical solution is given by Sedov 11 , von Neumann 12 , and 
later by J.L. Taylor 22 . Sedov 11 obtained solutions for plane, cylindrical and spherical blast waves, and for various 
values of y . Von Neumann 12 sought a similarity solution for the point explosion (treated by Taylor) in Lagrangian 
formulation. The unknown shock boundary as a free surface is determined as part of the solution. The solution for 
the blast wave problem in Lagrangian coordinates is explicit in terms of the parameter 9 = £ c / £ i (ratio of kinetic 
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and internal energies). This is in contrast to the implicit solution of the Eulerian equation of motion found by 
Sedov 11 and J.L. Taylor 12 , see (Sachdev 2000). 

Fig. 4a shows the prediction for the variation of the shock radius as a function of time, which is in close 
agreement with the measurements. The solutions (Sedov 11,8 ) are shown in Fig. 4b, and exhibit three main features. 
The velocity curve (j> rapidly becomes a straight line passing through the origin; 2. The density curve y/ 

approaches almost zero at t) « 0.5 , and remains close to it up to the center; 3. The pressure function f x decreases to 
become constant, and asymptotes to 0.37 times the maximum pressure just behind the shock. The very rapid decay 
of the density into the sphere is evident: almost all the gas is in a relatively thin layer behind the shock wave. This is 
due to the fact that the peak density of the gas immediately behind the shock ( 7] = 1 ) has density six times the 
normal density, as given by the Rankine-Hugonoit relations (Landau & Lifshitz 23 ). Most of the material in the blast 
wave gets accumulated near the shock. This tendency becomes more pronounced as y —> 1 . 

An approximate analytical solution, which is explicit and instructive, was reported by Taylor 5 introduced a 
correction to the linear behavior of u by adding a nonlinear term. Taylor 5 integrated the ODE in a closed form for 
the case of spherical blast wave. The resulting solution is much closer to the numerical solution. 

2. Shocks of Moderate Strength 

For intermediate times when the shock strength is finite, small departures from the classical solution due to 
counterpressure effects are accounted for in the perturbation solution by Sakurai 14 , and in the quasi-similar solution 
of Oshima 24 . As the shock strength decreases with time, there is seen a negative pressure phase (expansion region) 
behind the shock (Fig. 5a). Exact numerical solutions (Sedov 8 ) indicate an interesting feature, which is the 
development of a negative phase in the velocity distribution as the shock becomes weaker. Because the blast wave 

intensity decreases with time, the range ( X t < X < Xj ) of the blast wave profile increases, and the magnitude of the 

pressure peaks (maximum and minimum) of the blast wave decreases (Chen & Liang 25 ). Fig. 5b shows a typical 
pressure-time curve for the blast wave profile with a negative pressure region. 

The shock radius in the limiting cases of strong shock limit ( p 2 > p x ) and of acoustic limit corresponding to 
short and large distances are characterized by 


Rcc t 2/ 5 strong shock 
R = C 0 t sound wave limit (at large distances) 


(14a) 

(14b) 


a. Sakurai’ s solution 14,15,26 

Sakurai 14 generalized Taylor’s similarity solution, and proposed a perturbation scheme for moderately strong 
shocks, in which the Taylor solution for strong shocks appears as the zeroth order solution. Sakurai expressed the 
solution in terms of two variables (also see Sedov 8 , p. 297) 

4 = r/R(t\ V = h /17(f) f = cl IR 2 =\IM] (15) 


The Rankine-Hugonoit equations for shocks of arbitrary strength (eq. 5) are utilized here. As in the strong shock 
case, the energy released E a is assumed constant, and is expressed by 


E a = 


1 2 1 

—u + 

2 r - 1 
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P Po 


- pr a dr, a = 0,1,2 


(16) 
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where the a = 0 , 1 and 2 respectively corresponds to the plane, cylindrical and spherical blast waves. This 
represents the explosion energy per unit area of the surface of the shock when R equals unity. The Lagrangian form 
of the conservation of mass is 
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R a+i 


a + 1 


Thus the equation for the total energy becomes 
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(17a) 


(17b) 


The boundary conditions include the condition that the particle velocity at the center of explosion must be zero. 

Numerical and approximate solutions are reported by Sakurai 14 for the zeroth and first order solution. The 
accuracy of the approximate solution for a = 0 (plane) and 1 (cylindrical) for the zeroth order solution is 
remarkably good, the error being less than 5 percent (Sachdev 27 , p. 65). With regard to the first order solution, the 
expansion of Sakurai 14 provided inaccurate results for R —> oo that are not consistent with the exact numerical 
solutions (Sachdev 27 , p. 64). 

Sakurai 26 developed an approximate solution in a manner similar to that of Taylor 5,6 . The assumptions of 
Sakurai 26 are that the particle velocity profile behind the shock is linear and that the derivative of the density with 
respect to the shock Mach number in the continuity equation is small and may be ignored. This assumption 
introduces serious error in the weak shock regime, where the shock wave gradually decays to become a sound wave 
(Sachdev 27 , p. 64). The solution yields an asymptotic behavior for R —> oo, which is at variance with the well- 
known results obtained in this limit by Whitham 28 and Landau 29 . It is not in good agreement with numerical 

solution, covering the entire range 0 < (f> < 1 . For the asymptotic motion at large times (1/M 2 — > 1). Solutions 
have been obtained by Whitham 28 and Sedov 8 . 

The similarity parameter of Sakurai 9 is expressed by 

c 0 t/R 0 =f{R/R 0 ) (18) 


where C 0 denoted the undisturbed speed of sound, R the shock radius, and R 0 the characteristic length. A quasi- 
similarity rule was proposed by Oshima 24 . 

b. Solution of Bach and Lee (1970) 

Bach and Lee 17 reported an approximate solution similar to that of Sakurai 26 for the entire range of the blast 
evolution (from strong shock to weak acoustic wave). The main assumption is that the density behind the shock 
wave in the blast obeys a power law, with the exponent being a function of time. A quasi-similarity exists for the 
decay coefficient 0 for moderate blast strengths as (Higashino et al. 30 ) 

0 = O{rj) (19) 

where 0 = RR/R 2 =— 1^-1 T) = \/M ] (20) 

M s \ dR ) 
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As shown in Fig. 6a, the results for a spherical wave with y — 1 .4 for the shock decay confident 6 (defined by 
eq. 20) agree closely with the numerical solution of Goldstine and Von Neumann 16 . The condition 
0 — — 1 corresponds to the strong shock limit, and 9 — 0 to the acoustic (weak shock) limit. Also shown in this plot 
are Sakurai’s perturbation solution (Sakurai 14 ’ 15 , and the approximate linear velocity profile solution of Sakurai 7 . 

The predicted pressure and particle velocity profiles are presented in Figs. 6b and 6c respectively for various 
values of the shock strength. The velocity profiles suggest that the particle velocities can be negative in regions 
behind the blast wave in the moderate and weak shock regions (i.e., particles are directed towards the center of 
symmetry). These results are consistent with the exact numerical solutions of Sedov 8 . 

3. The Shock Wave at Large Distance 

The similarity solutions were related to the description of strong shocks in the early stage of development 
(Sachdev). In the other extreme- the nearly linear regime, when the shock has become rather weak, Whitham 28 
proposed a general theory to cover such flows. It describes attenuation of spherical shocks at large distances from 
the origin. The only assumption is that since the entropy changes are of third order in strength, the flow with weak 

shocks may be considered isentropic. Thus R = 7?(c 0 ,£) at large distances form the origin. A fraction of explosion 
energy which is degraded as heat is thus not available for doing work as the shock propagates. 

C. High Pressure Gas Explosions (Bursting Sphere) 

A more realistic problem than that of point explosion is the case of a unit sphere of high pressure gas at time 
t = 0 suddenly expanding into a homogeneous atmosphere (air) at rest. This problem was treated numerically by 
Brode 18,19 and analytically (approximately) by Friedman 31 . At sufficiently high pressure, an important phenomenon 
arises here due to the occurrence of a secondary shock that does not appear in the one-dimensional shock tube 
problem (Sachdev 27 ). 

When a highly pressurized sphere is suddenly released, an inward rarefaction wave in the high-pressure region 
and contact discontinuity appears between the rarefaction wave and the shock wave (Fig. 7). At the contact 
discontinuity, the pressures and velocities on the two sides are the same, but the density and temperatures on the two 
sides are different. The contact discontinuity initially moves out behind the main shock front but with a decreasing 
velocity. A second shock wave can arise between the contact discontinuity and the rarefaction wave (Liu et al. 32 ). In 
the 1-D shock tube problem, the main shock and the expansion wave come into instantaneous equilibrium being 
separated by a region of uniform pressure and velocity (Friedman 31 ). Physically, the reason for the formation of the 
second shock is that the high pressure gas must expand to lower pressures than those reached through an equivalent 
one-dimensional expansion, due to the expansion of volume. On account of this overexpansion, the pressures at the 
tail of the rarefaction wave are lower than those transmitted back by the main shock, and a compression, or a second 
shock, must be inserted to connect the two phases. Because the second shock wave occurs in the expansion region, it 
is expected to be rather weak initially and propagate outwards with the expanding gas. However, its strength 
increases with time and reaches a fairly high intensity in a short time. Soon, the second shock wave stops 
propagating outwards and comes to a halt before reverting to implode on the origin. 

A typical blast wave pressure profile is shown here (Fig 8). The major difference between a planar shock wave 
and a planar blast wave is that the flow properties behind the shock wave front are constant. A distinct characteristic 
of a spherical blast wave is a moving shock wave (jump in pressure) immediately accompanied by expansion waves 
behind the blast wave front, there is a blast wave profile region in which the flow properties are rapidly decreased 
due to expansion waves (Chen & Liang 25 ). Due to its 3-D expansion, the intensity of the blast wave always 
decreases with time and distance during its propagation. The overpressure phase is due to the shock front and the 
under-pressure phase is due to the expanding waves following the shock front. Both phases result in two visible 
shock peaks on the pressure-time plot. By contrast in shock tube experiments, the pressure downstream of a shock 
wave is often constant between itself and the following contact surface. 

1 . Brode 9 s Numerical Solutions 

Brode 18 simulated spherical blast waves, and considered two cases of practical importance: 1. strong shock, 
point-source explosions, 2. the sudden release of high pressure isothermal (both cold and hot) gas from spherical 
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enclosures. For the point explosion problem, he considered the exact solution of von Neumann 12 as an initial 
condition. Ideal gas equation of state is considered. Artificial viscosity, introduced by von Neumann and 
Richtmyer 33 is considered. 

Brode 19 , in a subsequent study, considered a blast wave from a spherical charge of TNT, with an equation of 
state for trichlorotoluene (TNT). In his numerical studies of detonation with a spherical charge of TNT, the existence 
of a second shock in the focal region was demonstrated. By considering real gas at high temperature, Brode 19 has 
shown that the blast over pressure at any radius is lower than that given by the ideal gas solution, since a part of the 
available energy is absorbed by the processes of ionization and dissociation. 

Goldstine and von Neumann 16 and Brode 18 show that the variation of shock overpressure varies with the shock 
radius R follows a power law of the form 


p-\ = AR n (21) 

where n is a slowly varying function of R , and A a constant (Sachdev 27 , p. 242). 
a. Brode ’s Solution for Blast Overpressure 

Brode 19 has shown theoretically that the peak overpressure p s0 from a free spherical air blast can be expressed 
in dimensionless form (Wyle Labs 34 ) 


P so r d 

Po {E/Po) 



=/W 


( 22 ) 


where p so is the peak overpressure due to TNT explosion, p Q the atmospheric pressure, X the reduced distance, 
d the actual distance (ft), and E the energy of the blast (lbf.ft). 

For standard atmospheric conditions, 


E_ 

Po 


PF r xl.55xl0 6 

2U6 


= 733 W T 


(23) 


where W T is the TNT equivalent weight (lb). Thus 

Pso=fz( Z ) (24a) 

where Z = d! = 9.02 X (24b) 

For surface explosions, the blast parameters are equivalent to a free air blast with twice the energy of free air 
blast. Thus 


Z = 2 1/3 9.02/1 = 11.37/1 

Fig. 9 represents a curve fit to Brode’s numerical solution for the blast overpressure from TNT explosions as a 
function ofZ . A power-law fit appears to be satisfactory. The numerical solution agrees with the test data on TNT 
and propellant explosions. 
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Fig. 10a shows the variation of shock Mach number and shock speed as a function of overpressure, as given by 
eq. (8a). The shock speed increases with increasing overpressure. For example, at an overpressure of 400 psia, the 
shock Mach number is about 5, with the shock speed being 5500 ft/s, which is comparatively high relative to the 
ambient sound speed (1100 ft/s). The dependence of particle velocity on overpressure is provided in Fig. 10b, as 
obtained form eq. (8b). 

D. Spectral Characteristics and Detection of Blast Waves 

Blast waves could be detected by a wideband microphone (Diaci 4 ), who investigated the weak shock region. By 
analyzing blast wave signals in the time and frequency domain , the characteristic features of their nonlinear 
waveform evolution can be investigated. Many blast wave studies were concerned with the kinematic analysis of 
shock wave propagation, including shock position vs. time or shock velocity vs. time. The decay of blast wave 
amplitude (pressure, density or any other physical quantity) has been studied less often. Typical detected signals and 
their amplitude spectra are shown in Fig. 11. The basic signal waveform shows an initial positive part that 
corresponds to compression followed by a negative part attributed to rarefaction. The spectra exhibit typical 
bandpass characteristics with one central peak and approximately constant percentage bandwidth. The higher the 

amplitude, the larger the signal duration t w . The higher the spectral amplitude peak, the lower the peak frequency. 


The strong dependence of the waveform duration on the blast energy E h is also reflected in the frequency 

domain. The spectral peak f m increases %sE h decreases. The nonlinear broadening of the waveform due to shock 

wave propagation is observed. As the spherical blast wave propagates, the tail of the rarefaction gradually steepens 
and eventually a shock builds up there also. Asymptotically, the waveform evolves into the well-known symmetrical 
N-waveform with shock discontinuities iii the front and the back. As a result of this distortion, ripples are introduced 
into the high frequency part of the spectrum. 

The following section outlines general (nonformatting) guidelines to follow, drawn from the original AIAA 
Manuscript Preparation Kit. These guidelines are applicable to all authors (except as noted), and include information 
on the policies and practices relevant to the publication of your manuscript. 

ID. Interaction of Spherical Blast Waves (Blast Wave Coalescence) 

Reported studies on interacting spherical blast waves are relatively few, although The interaction of a pair of 
(weak) cylindrical blast waves of either equal or unequal strength were studied experimentally and computationally 
by Higashinoet al. 30 . In sequential detonation and spatially distributed blast sources, it is known that shock waves 
coalesce and focusing can lead to significant amplification of peak overpressure. At sufficiently large distances from 
the charges, coalescence will occur. An overtaking wave may first encounter the negative overpressure phase of the 
wave being overtaken, causing it to decelerate. The overtaking wave must be somewhat stronger than the wave 
being overtaken or it will not be able to penetrate the negative overpressure phase and coalescence cannot occur. 
Once this phase has been penetrated and the overtaking wave has penetrated the positive phase of the wave being 
overtaken, coalescence takes place but occurs only after some additional propagation. 

Several approximate methods of combining peak overpressure from N coalesced waves are described in 
Starkenberg and Benjamin 1 as follows: 

A. Simple superposition 


i Y 


( 22 ) 


where A p represents the combined overpressure, and tSp x the overpressure of blast wave i . 

B. Full LAMB Model 


Hikida and Needham 35 utilized the following scheme for blast coalescence. 
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P = A+Z a A- 

i=l 

1/ =-!=!— (23) 

1=1 ^ v z /=i ^ / 

where p refers to density, and V the velocity of blast propagation. 

C. Single Effective Charge 

The method of single effective charge considers the use of pressure produced by the total explosive weight of 
all stacks with the coalesced waves combined at their centers of charge. That is, the peak overpressures associated 
with coalesced waves are essentially the same as those produced by single charges of the same total explosive 
weight. Zaker 36 experimentally showed that with three equal charges, the third pulse tends to overtake the second 
before the second overtakes the first. 

Charges of equal/unequal mass (charges) and delay time, and separation distances need to be considered in the 
development of pressure combination algorithms and coalescence maps (patterns). 

IV. Test Data on Overpressure 

Blast overpressure data were obtained using scale model and full scale solid rocket propellant motors (upper stage) 
at NASA White Sands Test facility 37 . The solid propellant is struck by shaped charges of 500 gm and 20 gm. Both 
symmetrical (equal strength) and unsymmetrical (unequal strength) wave interactions were studied. 

A. Scale Mode! Tests 

1. Test Configuration and Instrumentation 

In the scale-model test, considered here, the propellant weight was 474 lb. For comparison, the full scale 
propellant in the full scale motor is about 4,431 lb. The configuration (Fig. 12a) includes a Ti-cylinder around the 
propellant stack to simulate the effects of the casing around the actual motor. The solid propellant was impacted 
with one 500-gm conical-shaped charge of C4 and 17 gm of Comp A-4 booster pellet. The burning of notable 
unbumed propellant fragments that were collected is videotaped and timed. The data is recorded up to 45 sec. from 
the fire command (t = 0 sec). 

The blast .overpressure is measured by piezoelectric transducer pressure gauges (PCB) mounted on three 
adjustable elevated gauges lines (EGL), namely F, G and H, separated 120 deg. apart. EGL F and H are fitted with a 
flat plate at one end, with the PCB mounted in the plate center (thus measuring the stagnation overpressure, which 
includes both the static and dynamic components of overpressure (Fig. 12b). Five pairs of PCB’s are flush mounted 
equal distances apart on each side of the pipe, in designated positions A and B. The PCB’s are located 180 deg from 
each other, parallel to the ground. The EGL’s are constructed from 4 in. diameter steel pipe, 30 ft. long, and are held 
72 in. above the ground (parallel to). The PCB’s have a 500 kHz resonant frequency, and a 1-//S rise time. Only the 
near-field blast pressure is measured. 

2. Test Data for Overpressure 

Fig. 13a shows the overpressure history for the stagnation overpressure from the closest gauge on EGL-F 
(channel-67), which is 4.92 ft from ground zero. This typical blast waveform is characterized by a peak overpressure 
(30 psia) at 16.2 ms, and positive pressure duration (about 1 ms). The stagnation overpressure includes a dynamic 
component due to the particle motion. We also notice a distinct wave due to propellant explosions, with a peak 
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overpressure of 16.4 psi at 19.2 ms. Only one major peak overpressure due to the propellant bum is observed here. 
The peak overpressure due to the charge is higher than that due to that of propellant. Both the blast overpressure and 
the positive pressure duration are consistent with the blast scaling due to Brode (1950). It is remarkable that the 
propellant pulse is wider than the pulse resulting from the 500-gm CSC. The relatively longer decay of the 
overpressure from the propellant blast is attributable to the long duration over which the energy is added to the 
flowfield. The overpressure phase is associated with the shock front location, and the under-pressure phase is linked 
to the expansion waves following the front. 

The overpressure history for a number of locations along the EGL: F-A is demonstrated in Fig. 13b. This 
illustrates the propagation and attenuation of the blast wave. For example, considering the first two PCB’s 
(separated by about 2 ft), the 500-gm peak is displaced by about 1 ms. This result is consistent with the shock speed 
of about 1500 ft/s (Fig.lOa). The peak A p diminishes from 30 psi at 4.9 ft to 2.5 psi at 23 ft. 


Fig. 13c presents an extrapolation of the near-field Ap data from the near-field to the far- field covering 
distances (altitudes) from 25,000 ft. The Ap decay seems to be somewhat consistent with that indicated by the 
scaling laws. 

B. Full Scale Propellant Test 

1. Test Configuration and Instrumentation 

The full scale tests are described in WSTF Special Test Report 37 . Fig. 14a shows a photographic view of the 
actual motor. The overall dimensions of the test article were 5 ft wide x 7.5 ft tall with 4431 lb of propellant. The 
propellant composition is: 71 % ammonium perchlorate, 18 % aluminum powder, 11 % hydroxyl-termicrated 
polybutadiene binder (HTPB). The motor casing is fabricated from titanium alloy TI -6A1-4V. A total of four 
shaped charges (CSC) were used. An aft enclosure was used to simulate the nozzle. Two 500-gm shaped charges 
were used for the breakup system, mounted in the PAF oriented 180° from each other. Two 20-gm shaped charges 
were located in the spin table assembly, near the top of the nozzle skirt. The two 20-gm shaped charges were 
oriented 180° from each other in the same vertical plane as the 500-gm charges. For the PNH BUS, the 20-gm 
charges were located at azimuths of 191.3° and 338°, while the 500-gm charges were at 90° and 270°. All shaped 
charges were initiated in the timing sequence designed for the MER destruct system. In the MER BUS design, the 
500-gm CSCs are initiated 0.87 ms after the 20-gm shaped charges are actuated. The PNH BUS timing sequence is 
similar. 

The test set up and instrumentation is essentially similar to those described for the scaled propellant tests (Fig. 
14b). The instrumentation for the full scale testing is similar to that in the scaled tests, except the charges are not 
symmetric for EGL-F and EGL-G. 

2. Test Data for Overpressure 

Fig. 15a illustrates the history from the stagnation transducer on EGL-F (channel-67). It indicates a series of 
blast waves. Two 20-gm charges are followed by two 500-gm charges in succession. The peak Ap due to 500-gm is 
considerably higher than that due to 20-gm peaks, since the energies are different, as indicated by the scaling law 
(Appendix-A.l). For the EGL-F the two 500-gm are separated by about 4 ft., and travel at about 2300 ft/s, so that 
At « 1.7 ms (the measured At is about 1.7 ms). The peak Ap related to propellant explosion are superposed on to 
the 500-gm charge waveforms, and the peak magnitude (about 20 psi) are consistent with those for the scaled model 
test. Also secondary shock waves are noted. The 20-gm charges do not seem to have initiated any significant 
propellant bum (or breakup), with more localized rather than large-scale action. 

The peak Ap history for various locations along the EGL-F (Fig. 15b) suggests that the 500-gm charges 

ultimately coalesce at about 20 sec. This merging of shocks is a consequence of the fact that the trailing shock 
travels relative to the particle motion occasioned in the wake of the leading shock; thus it moves at a higher speed 
for the same strength. The Ap decreases from 65 psi to 7.5 psi as the distance increases (due to spherical spreading). 

The strength of the merged shock is due to the two charges. There seems to be a problem (data drop out) with 
channel 37 after about 12 sec. 
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The characteristics of the blast propagation for EGL: F-B (B-side transducer) is seen to be similar to that of 
EGL: F-A in all respects (Fig. 15c). This is to be expected, as the transducers are similarly situated. This 
demonstrates the consistency of measurements. Likewise, in the case of EGL-G (Fig. 15d), the behavior is similar to 
that of EGL-F, since the charges are similarly situated. In the case of EGL-H (Fig. 15e), the two 500-gm charges are 
located equidistantly relative to the gauges, so that the two peaks are very close. 

The spectral distribution of energy for the pressure signal from channel 67 (F-A) for individual explosions of 
short duration for the various explosions (20 gm charge, 500 gm charge, and propellant explosion) are depicted in 
Fig. 16. It is seen that the energy is distributed over a broad range of frequencies, as is to be expected for impulse 
type functions typical of the blast wave. There is some spectral similarity for the various explosions, with the spectra 
displaying a distinct peak. The spectra are similar to those presented for blast waveforms (Diaci and Mozina 1992). 

Fig. 17 presented the extrapolation of peak pressures to the far-field. The extrapolation suggests that the peak 
pressures at about 25,000 ft are discernible, thus pointing to the feasibility of the ground microphone system for on- 
board explosion detection. However, the occurrence of multiple explosions, blast wave coalescence, and absorption 
of sound at high frequency and over long distances should be taken into account. 


V. Conclusions 

Both the scaled test data and full scale propellant blast test data have been helpful in our understanding of the blast 
propagation from single and multiple charges. The data illustrate the effects of blast coalescence and propagation of 
blast waves of different strength. An examination of the test data suggests that the test data for blast overpressure are 
generally consistent with the available blast scaling laws. Further work is in progress with regard to the interaction 
effects of blast waves. 
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Fig. 1 Temporal evolution of a blast wave 



Fig. 2a Schematic of a point source explosion 
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Fig. 2b Schematic of an explosion of an initially spherical charge. 
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Fig. 3 Propagating and stationary shocks (adapted from Liepmann & Roshko, 1957). 
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Fig. 4a Variation of shock radius with time for the strong shock solution of Taylor (1946). 



Fig. 4b The strong shock solution of Sedov (1946). 
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Fig. 5a Spatial distribution of pressure at various times for intermediate shock strength 
(adpated from Kinney 1962). 



Fig. 5b Pressure-time varaition of a blast wave of moderate shock strength 
(adapted from Kinney 1962). 
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Fig. 6a Comparison of shock decay coefficient for shocks of moderate strength 
(Bach & Lee 1970). 



Fig. 6b Pressure distribution behind the shock for various times (Bach & Lee 1970). 
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Fig. 6c Particle velocity distribution at various time (Bach & Lee 1970). 
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Fig. 7 Comparison of shock positions for the explosion of a spherical charge (Friedman 
1971). 
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Fig. 8 A typical blast wave from a spherical charge, showing secondary shock. 



Fig. 9 Brode’s numerical solution for shock overpressure (Brode 1959). 
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Fig. 10a Shock velocity and shock Mach number from Rankin-Hugonoit reations. 
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Fig. 10b Particle velocity for blast wave propagation. 
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Fig. 1 1 Spectral characteristics of blast waves (Diaci 1992). 
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Fig. 12a Test set up for scaled testing. 



Fig. 12b EGL and PCB locations for scaled test. 
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Fig. 13a Overpressure history at the stagnation location for the scaled test. 
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Fig. 13b Overpressure history for the scaled test (contd.). 
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Fig. 13c Peak overpressure scaling for the scaled test data. 
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Fig. 15a Overpressure history at the stagnation point for full scale testing. 
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Fig. 15b Overpressure history of sensors along F-A for full scale testing. 
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Fig. 15c Overpressure history of sensors along F-Bfor full scale testing. 
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Fig. 15d Overpressure history of sensors along G-A for full scale testing. 
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Fig. 15e Overpressure history of sensors along H-A for full scale testing. 



Fig. 16 Spectral distribution for full scale testing. 
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Fig. 17 Peak overpressure scaling for full scale testing. 
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